Mon. Not. R. Astron. Soc. 000, [TV?? (2005) Printed 5 May 2009 (MN LM^X style file v2.2) 



RS Ophiuchi: Thermonuclear Explosion or Disc Instability? 
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ABSTRACT 

Sokoloski et al (2008) have recently reported evidence that the recurrent nova RS Ophi- 
uchi produced a pair of highly collimated radio jets within days of its 2006 outburst. 
This suggests that an accretion disc must be present during the outburst. However in 
the standard picture of recurrent novae as thermonuclear events, any such disc must 
be expelled from the white dwarf vicinity, as the nuclear energy yield greatly exceeds 
its binding energy. We suggest instead that the outbursts of RS Oph are thermal- 
viscous instabilities in a disc irradiated by the central accreting white dwarf. The 
distinctive feature of RS Oph is the very large size of its accretion disc. Given this, it 
fits naturally into a consistent picture of systems with unstable accretion discs. This 
picture explains the presence and speed of the jets, the brightness and duration of the 
outburst, and its rise time and linear decay, as well as the faintness of the quiescence. 
By contrast, the hitherto standard picture of recurrent thermonuclear explosions has 
a number of severe difficulties. These include the presence of jets, the faintness of 
quiescence, and the fact the the accretion disc must be unstable unless it is far smaller 
than any reasonable estimate. 
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1 INTRODUCTION 

Sokoloski et al (2008) have presented evidence that within 
days of its 2006 outburst the recurrent nova RS Ophiuchi 
produced a pair of highly collimated ('half-opening angles 
of a just few degrees') and high velocity (v ~ 5000 km s -1 ) 
jets. Taylor et al. (1989) made similar deductions after the 
previous (1985) outburst. We argue here that this provides 
compelling evidence about the cause of the outburst in this 
object. 

The source of the energy driving the outbursts in re- 
current novae has been under discussion for some time (e.g. 
Livio, Truran & Webbink, 1986). There are two basic possi- 
bilities: either the outburst is caused by a nuclear explosion 
on the surface of a white dwarf, akin to a classical nova ex- 
plosion, but smaller in scale, or the outburst is powered by 
accretion energy released by a change in the accretion rate 
on to the central white dwarf, akin to dwarf nova outbursts 
but larger in scale, and similar in nature to the transient X- 
ray binaries (which differ in having a central neutron star or 
black hole as the accretor). 

It is, however, difficult to see how a nuclear-powered ex- 
plosion on the surface of an accreting white dwarf could give 
rise to such well-collimated jets. In general one would ex- 
pect nuclear burning in the hydrogen-rich surface shell to be 
quasi-spherically symmetric. The nuclear energy yield from 



the hydrogen-rich matter considerably exceeds its binding 
energy at the surface of a white dwarf. In a nuclear-powered 
nova, the mass of explosively burning material is compara- 
ble to that accreted since the last outburst. The mass in the 
central regions of any disc must be much smaller than this, 
and is even less gravitationally bound than the matter on 
the white dwarf surface. Thus the mass in such a disc does 
not have enough inertia to provide strong collimation of the 
flow. This point is underlined by the lack of any strongly 
collimated flows to be seen in ejecta from classical novae 
(Slavin et al., 1995; Gill & O'Brien, 1998). One only sees 
slight asymmetries in the shells of material expelled by the 
nova explosion, presumably caused by the presence of a close 
binary companion, and lumpiness in the shells (presumably 
caused by instabilities in the ejection process). 

In contrast, the jets seen in RS Oph are high velocity 
and strongly collimated. As such they are strongly reminis- 
cent of jets seen in other astronomical objects such as radio 
galaxies, active galactic nuclei (AGN), microquasars, some 
binary X-ray sources and young stellar objects. These jets 
all have the following in common (Pringle, 1993; Livio, 1999, 
2000; Price, Pringle & King, 2003): the source of the jet has 
accretion occurring through a disc; the jet velocities are com- 
parable to, or slightly higher than, the escape velocity from 
the central accreting object; and about 10 per cent of the 
accreting material is being ejected in the jet. There is also 
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a suggestion that the collimation of such jets require a large 
ratio of inner to outer disc radii. It is immediately clear that 
the jets seen in RS Oph fall precisely into this category. 

We therefore argue here that we should abandon the 
view that the recurrent outbursts of RS Oph are nuclear- 
powered. We propose instead that the outbursts are driven 
by gravitational accretion energy, so that the system is ac- 
tually a dwarf nova with unusual properties which all stem 
from the fact that its accretion disc is far larger in size than 
in normal dwarf novae. This property leads to unusually 
bright, long-lasting outbursts, in which irradiation of the 
disc by the central source probably plays a role, and simi- 
larly long quiescences. Many such outbursts are observed to 
have jets. 

In the rest of this paper we first show that a dwarf nova 
origin offers a reasonable explanation of the outbursts of RS 
Oph, and then examine the evidence for a thermonuclear 
origin. Section 4 is a discussion. 



2 RS OPH AS AN IRRADIATED DWARF 
NOVA 

The orbital period of RS Oph is 453.5 days (Brandi et 
al., 2009; Fekel et al., 2000), implying a binary separation 
a ~ 1.4 AU. The companion star is a red giant, and the 
measured period gives its orbital velocity as 30Mj L//3 (l + 
M 2 /Mi)- 2/3 kms" 1 where M 1 ,M 2 are the white dwarf and 
giant masses in Mq, and mass transfer stability requires 
M2 < Mi . (Assuming that the emission wings of the Ha line 
trace the motion of the white dwarf, Brandi et al (2009) 
give Mi — 0.59 sin -3 i, M2 = 0.35 sin -3 i, with i the inclina- 
tion.) Mass transfer could potentially occur either through 
capture of the red giant wind or Roche lobe overflow. How- 
ever since the stellar wind velocity ~ 20 kms -1 (Bode & 
Kalm, 1985, p. 206) is unlikely to exceed the giant's orbital 
velocity the dynamics of the mass transfer must be simi- 
lar to Roche lobe overflow even if the giant's photosphere 
does not fill this lobe. The system thus has an accretion disc 
with outer radius ~ 70% of the white dwarf Roche lobe, 
hence size ~ 0.5a ~ few xlO 12 cm. Unless the outer radius 
of the accretion disc is far smaller than this estimate, i.e. 
~4x 10 10 cm ~ 2x 10~ 3 a, the gas there must be much 
cooler than the typical temperature Th ~ 6500 K for hy- 
drogen ionization, implying that the system must undergo 
disc instabilities. 

A disc instability at radius R ~ 10 12 i?i2 cm has charac- 
teristic viscous timescale t v isc ~ R 2 /v with v the kinematic 
vicosity. We adopt the standard parametrization v = ac s H , 
where c s is the local sound speed and H the disc scale- 
height, with a a dimensionless parameter. In a thin, fully- 
ionized disc, a wide range of observations give a typical range 
a ~ 0.1 - 0.4 (King et al., 2007). Accordingly for the hot, 
outburst state of the disc in RS Oph we take the viscosity 
coefficient as qj, = O.lao.i. Using eqs (5.49) of Frank et al., 
2002, this gives 

tvisc = lSOo^^R^M^ 3 days (1) 

where R12 is the disc radius in units of 10 12 cm, and M21 
the peak outburst accretion rate in units of 10 21 gs -1 . Im- 
mediately before the outburst begins, a significant part of 
the disc must reach the maximum surface density E max (i?) 



for a range of disc radii R, and one could in principle es- 
timate the unstable disc mass by integrating this quantity. 
But E m ax(-R) depends on the cooZ-state disc viscosity, which 
is essentially unknown, so we take instead E(i?) = /E m j n = 
8.25 x 10 2 /i?i2 05 M 1 -°- 35 a^ - 8 gem" 2 . Here E min is the min- 
imum surface density attained at the end of the outburst, 
so that / > 1 (cf eqn (5.84) of Frank et al., 2002). Integrat- 
ing this over the disc surface we can estimate the total disc 
mass Mh before outburst. We approximate the local peak 
accretion rate as M ~ 2Af/ l /t v i sc , as appropriate for an ex- 
ponential or linear viscous decay (see below) and using 
above we can now solve for M and t v i sc as functions of R12, 
i.e. 

t visc = 127a^?- 8 R° 12 i8 M?- 15 days (2) 
and 

M pcak = 2 x lO^ipMr - 5 / 1 ' 43 gs -1 (3) 
with the consequence 

»> c c w in 2 0j.5.4 4.3,^-1/3^.1.43 -1 / A \ 

M peak = 5.5 x 10 t 10O a , 1 M 1 f gs (4) 

where tioo is tvisc in units of 100 days. 

Before using these relations for RS Oph, we note that 
they give very reasonable estimates for known cases of vis- 
cous decays if we take the parameter / = E max /E m i n ~ few, 
as indeed suggested by numerical modelling of disc instabil- 
ities (see Lasota, 2001) for a review). 

Thus for a superoutburst of an SU UMa system we 
take Mi = 0.7, R12 = 2.5 x 10~ 2 and find a peak out- 
burst accretion rate 2x 10 17 / 1 ' 43 gs -1 (accretion luminosity 
L ~ 10 34 / 1 ' 43 ergs -1 ) and a decay timescale of ~ 14 days. 

For a neutron-star X-ray transient we take R12 = 
0.1, Mi = 1.4 and find a slightly super-Eddington luminos- 
ity L ~ 5 x io 38 / 1 43 ergs" 1 with t viBC ~ 42 days. 

Both estimates, and others that one can make for e.g. 
black-hole accretors in transients, agree well with observed 
systems, so we now apply eqns ((2] |3j) to RS Oph. 

2.1 Application to RS Oph 

For RS Oph we take R12 — 1, Mi ~ 1 and find a peak 
accretion rate 

Af poak ~ 1.7 x 10 21 / 143 gs" 1 = 2.5 x lO^Meyr" 1 , (5) 
with an outburst timescale 

t v i sc — 128c%i' 8 days. (6) 

The accretion rate © is somewhat super-Eddington 
(for gravitational energy release), formally giving a lumi- 
nosity GMM pei: ,k/R = 4.6 x 10 38 ergs -1 for accretion on to 
a IMq white dwarf, so that M pe ak — 3A'/ E dd- In line with 
expectations from other super-Eddington accretors such as 
ultraluminous X-ray sources we anticipate strong outflow 
from the inner parts of the accretion disc, leading to an ac- 
cretion luminosity L ~ L E dd[l + ln(A'/p Cak /Af E dd)] = 2L Edd 
(cf Shakura & Sunyaev, 1973; Begelman et al., 2006, Pouta- 
nen et al., 2007). At this value (M pea k/MEdd — 3 the emis- 
sion is almost isotropic (at higher values of this ratio it is 
collimated towards the disc axis). The system should thus 
have a luminosity of ~ 2L E dd unless viewed at very high in- 
clination. This agrees with observations of RS Oph, as does 
the outburst duration implied by (|6}. 
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The rise of RS Oph to maximum was very rapid 
( < 1 day) . This follows from the dwarf nova picture given 
here also. To brighten from an accretion luminosity of 
10 37 ergs -1 to 10 38 ergs -1 requires the disc instability heat- 
ing front to move from disc radius Rd ~ 10 10 cm to the 
white dwarf radius, which is about 10 times smaller. Eqn 
(51) of Lasota (2001) shows that the disc instability must 
be of 'inside-out' type, as the opposite would require a mass 
transfer rate greater than about 4 x 10 20 gs -1 , i.e. compa- 
rable with the outburst accretion rate. For such outbursts 
the front moves with a velocity < cthC s , so the rise time is 
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(7) 



Using the disc central temperature given by Frank et al., 
2002 (eqn 5.49) to estimate c s we find t T i sc > 3 hours. 

Thus far the dwarf nova model gives a consistent repre- 
sentation of the outburst of RS Oph. However we need one 
modification of the usual picture. A normal dwarf nova out- 
burst would end not after a viscous time, but after a local 
thermal time, which is considerably shorter. The crucial dif- 
ference in RS Oph is the very large size of the accretion disc. 
This means that irradiation by the central source (temper- 
ature oc i? -1 / 2 ) must dominate local viscous energy release 
(temperature oc P -3 / 4 ) at large disc radii, where most of the 
disc mass is. The outburst then resembles those of soft X-ray 
transients: the disc is trapped in the hot outburst state un- 
til central accretion drops because most of the heated mass 
has been accreted (King & Ritter 1998). For a disc heated 
by a central point source of luminosity L c the irradiation 
temperature Ti rr (P) is given by (cf King & Ritter, 1998) 
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where a is the Stefan-Boltzmann constant, g is a geometric 
factor of order 0.2, and H/R ~ 0.04 is the disc aspect ratio. 
This gives an irradiation temperature of order 
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that is, the central accreting white dwarf keeps the disc self- 
consistently in the hot state (Ti rr > Th) out to a radius 
Rh ~ 10 12 cm at the start of the outburst. In this case the 
outburst must evolve on the local hot-state viscous time, 
and decay exponentially or linearly in time depending on 
whether irradiation keeps the whole disc in the hot state or 
not (King & Ritter, 1998). As Rh is smaller than the full 
disc outer radius, the outburst decay should be (bolometri- 
cally) linear over a viscous timescale. This is consistent with 
observation (Page et al., 2008, Fig.l) between days 50 to 
100. 

Between outbursts we expect the disc to refill, with 
little or no accretion on to the white dwarf. This offers a 
natural explanation for the faintness of the system between 
outbursts, in particular in X-rays (Mukai 2009). We con- 
clude that the main features of the outbursts of RS Oph 
are successfully explained in terms of an instability in a disc 
irradiated by the central accreting white dwarf. 



3 RS OPH AS A RECURRENT 
THERMONUCLEAR SOURCE 

We have suggested above that the outburst behaviour of 
RS Oph is explicable in terms of instabilities in an irra- 
diated disc. Until now the usual interpretation of RS Oph 
and other recurrent novae is that matter accumulates on 
the surface of the white dwarf until it reaches the base pres- 
sure Pcrit ~ 2x 10 19 dyne cm -2 required to ignite as a 
runaway thermonuclear event. We note that the presence of 
processed material in the outbursts is not per se proof of 
nuclear burning in the outburst, since the companion star is 
evolved. Thus any material transferred to the white dwarf 
and subsequently expelled in an accretion-driven outburst 
came originally from the convective envelope of the red gi- 
ant, where it mixed with products of the nuclear-burning 
shell around the degenerate core. Its composition therefore 
does not differ greatly from that of matter partially burnt 
and ejected from the white dwarf surface in a classical nova. 
This picture has several difficulties, which we list below. 

1. The presence of jets requires an accretion disc, and is 
very hard to reconcile with a thermonuclear model, as noted 
in the Introduction. The nuclear energy yield from burning 
hydrogen-rich matter considerably exceeds its gravitational 
binding energy at the surface of a white dwarf. Hence a ther- 
monuclear explosion would inevitably blow away the inner 
parts of any accretion disc. 

2. Even assuming a mass M\ close to the Chan- 
drasekhar limit, the white dwarf has to accrete at a rate 
M = Pcrit^/GMitqu > 1O -8 M yr -1 in order to accumu- 
late enough mass to ignite the outbursts (here Pi is the 
white dwarf radius, and t qu ~ 20 yr is the duration of quies- 
ence). But the accretion rate M>10 -8 Mq yr -1 is hard to 
reconcile with the lack of X-rays observed from RS Oph in 
quiescence. Mukai (2009) considers various possible ways out 
of this conclusion, including combinations of very high in- 
trinsic absorption columns and totally optically thick bound- 
ary layers, and concludes that none are convincing. In par- 
ticular, CVs with boundary layer emission always in practice 
have a surface layer optically thin enough to produce hard 
X-rays (Patterson & Raymond, 1985), i.e. an optically thin 
region where accretion energy is released. 

3. The very wide binary orbit of RS Oph suggests that 
its disc is likely to be unstable, and thus reduce central ac- 
cretion severely below the value M > 1O -8 M0 yr -1 required 
for the thermonuclear model. To avoid this the outer disc ra- 
dius must be smaller than ~ 2 x 10 -3 a, where a is the binary 
separation. This is ruled out if the giant fills its Roche lobe, 
and very unlikely even if mass transfer is via stellar wind 
capture, since the wind speed is rather less than the orbital 
velocity of the giant companion. 

We note finally that the dwarf nova-type outbursts con- 
sidered in Section 2 imply a mean white dwarf accretion rate 
Mo < 7 x 10 -8 Mq yr - 1 if they recur at the current rate of one 
per ~ 20 yr (note that mass loss during the outburst makes 
this an upper limit) . With a non-extreme white dwarf mass 
one might then speculate on a thermonuclear nova occur- 
ring after a few outbursts. However the ultimate source of 
mass is not the disc, but the companion star. Unless this 
transfers mass at a rate > Mo, the disc will run out of mass 
and the outbursts will recur more slowly. This is another 
way of saying that (as always) the nova recurrence time is 
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~ AM/ (— M2), where AM is the mass that has to accumu- 
late on the white dwarf surface. For mass transfer on the nu- 
clear timescale of a low-mass giant (as here), the likely rate 
—M2 is less than Mo by factors > 3 (cf Webbink, Rappaport 
& Savonije, 1983). This suggests a 'supercycle' of outbursts 
and a nova recurrence time > 10 3 yr. This is short compared 
with most CVs simply because — M2 is much larger here. 
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4 DISCUSSION 

The distinctive feature of RS Oph is the very large size of its 
accretion disc. Given this, it fits naturally into a consistent 
picture of systems with unstable accretion discs. This picture 
explains the presence and speed of the jets, the brightness 
and duration of the outburst, and its rise time and linear 
decay, as well as the faintness of the quiescence. We would 
expect this general picture to apply to other long-period 
systems, e.g. T CrB. 

By contrast, the hitherto standard picture of recurrent 
thermonuclear explosions has a number of severe difficulties. 
These include the presence of jets, the faintness of quies- 
cence, and the fact the the accretion disc must be unstable 
unless it is far smaller than any reasonable estimate. 

A simple observational test of the thermonuclear picture 
would be provided by a dynamical mass determination of the 
white dwarf mass. A thermonuclear origin for the outbursts 
requires a mass extremely close to the Chandrasekhar value 
~ 1.4M©. 
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